The Hamiltonian description of incompressible fluid ellipsoids 



P.J. Morrison* 

Department of Physics and Institute for Fusion Studies, The University of Texas at Austin, Austin, TX, 

78712, USA 

Norman R. Lebovitz 
Department of Mathematics, The University of Chicago, Chicago, IL, 60637, USA 

Joseph A. Biello 

Department of Mathematics, University of California at Davis, USA 



Abstract 

We construct the noncanonical Poisson bracket associated with the phase space of first 
order moments of the velocity field and quadratic moments of the density of a fluid with a free- 
boundary, constrained by the condition of incompressibility. Two methods are used to obtain 
the bracket, both based on Dirac's procedure for incorporating constraints. First, the Poisson 
bracket of moments of the unconstrained Euler equations is used to construct a Dirac bracket, 
with Casimir invariants corresponding to volume preservation and incompressibility. Second, 
the Dirac procedure is applied directly to the continuum, noncanonical Poisson bracket that 
describes the compressible Euler equations, and the moment reduction is applied to this 
bracket. When the Hamiltonian can be expressed exactly in terms of these moments, a 
closure is achieved and the resulting finite-dimensional Hamiltonian system provides exact 
solutions of Euler's equations. This is shown to be the case for the classical, incompressible 
Ricmann ellipsoids, which have velocities that vary linearly with position and have constant 
density within an ellipsoidal boundary. The incompressible, noncanonical Poisson bracket 
differs from its counterpart for the compressible case in that it is not of Lie- Poisson form. 
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1. Introduction 

The Euler equations governing tlie velocity field v, density p and pressure p of an inviscid 
fluid are 

^fV + v-Vv = -p"Vj9 + f and (1) 
aip + div(pv) = (2) 

where f denotes an as yet unspecified force. These equations must be augmented by boundary 
and initial data, and by further conditions relating the variables v,p and p: either an 
equation of state (in the compressible case) or the condition div v = (in the incompressible 
case). We shall be interested in exploring their Hamiltonian structure in a particular context. 
Our principal reference for a general discussion of this structure and the derivations of the 
corresponding brackets will be pQ. 

Exact solutions of the Euler equations are possible only under simplifying assumptions 
and in simple contexts. A family of solutions in the context of astrophysics, namely, where 
the force term f includes the self-gravitational effects of the fluid mass, exists under the 
assumption of a fluid of uniform density confined to an ellipsoidal domain, with a velocity 
field linear in the coordinates. These assumptions reduce the Euler equations to a finite 
system of ordinary differential equations. The equations for these Riemann ellipsoids have 
been widely investigated: their study goes back to the work of Dirichlet ([2]) and of Riemann 
([3]), but our principal reference for this will be [4J. We summarize their properties in 
Appendix |A| which we will refer to often. (In Appendix [B] we describe four natural frames 
of reference for the ellipsoids, which are included here because they do not appear to have 
been published together elsewhere.) 

In none of these references was the Hamiltonian nature of the finite-dimensional system 
emphasized. This was first addressed by Rosensteel [S] . His starting point was the so-called 
virial method originally introduced to investigate the stability of steady solutions of the 
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Euler equations. The virial is a moment of the form 



M-ij = / pXiVj d X = / XiMj d X, (3) 
Jd Jd 

where i,j = 1,2,3, and the second form introduces the specific momentum of the fluid 
M = pv. This moment is considered together with another moment, equivalent to the 
moment-of-inertia tensor, 

Ejj = / pXiXj d^x . (4) 
Jd 

Rosensteel presents an algebra for these moments, i.e., bracket relations among them that are 
closed, and that provide a noncanonical Hamiltonian description of the Riemann ellipsoids 
with a certain choice of the Hamiltonian function H^E, Ai); we present these relations below, 



in equations (12)-(14). 



We call attention to two features of Rosensteel's description of the incompressible case: 

1. The bracket relations are presented without reference to the fluid-dynamics equations 
g and (|2| ab ove, and 

2. The formulation requires a Hamiltonian function other than the total energy as well as 
the imposition of extraneous constraints. 

The feature ([T]) is addressed in ^ below, where we derive Rosensteel's bracket relations 
in a straigtforward way via a moment reduction of the general fluid-dynamical bracket ([T]). 
Feature ^ is discussed in detail in §[|5]and|6} 

We view the fluid as incompressible. This is natural because the density of the Riemann 
ellipsoids is spatially uniform]^ However, Rosensteel's bracket does not constrain the fluid 
to be incompressible, and we therefore modify it via Dirac's procedure for incorporating 
constraints. Dirac's method is described in ^ We observe in ^that one can alternatively 



^ There are also applications allowing for compressibility wherein p — p(t), i.e., the density is spatially 
uniform but varies with time. We do not address these cases here. 
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first apply Dirac's procedure and subsequently effect a moment reduction, with the same 
result. The resulting Dirac bracket is no longer of Lie-Poisson type: the bracket relations 
depend nonlinearly on the moments. 

In ^ we relate the noncanonical Hamiltonian equations obtained from Rosensteel's 
bracket to the equations describing the Riemann ellipsoids and show that, if the Hamil- 
tonian is taken to be the total energy, the pressure term from fiuid dynamics is missing. It 
can be restored by adding an extra term to the Hamiltonian. In ^we show that the Hamil- 
tonian equations obtained from the Dirac bracket using the total energy as Hamiltonian give 
the full equations for the Riemann ellipsoids, and, moreover, avoid the necessity of imposing 
any further constraints. Finally, in ^we summarize and discuss these results. 

2. The Lie-Poisson Bracket and its Moment Reduction 

The Euler equations ([T]) and (|2| can be re-expressed in terms of the momentum density, 
M := pv, as 

dtM + V ■ VM + (div v) M = - Vp + pf and (5) 
dtp + diYM = 0. (6) 

These, like equations ([T]) and ([2]), will be referred to as unconstrained, since neither the 
constraint of incompressibility nor that of an equation of state has yet been imposed. 

The Hamiltonian description of these equations is reviewed in The noncanonical 
Poisson bracket, as given in f6], i^ 

f f6G ^6F 6F „6G\ ,3 
+ / P t^-V — - — - V— d^x. 7 



■^Unless otherwise indicated, repeated Latin indices are summed from 1 to 3. 



4 



This is a Lie-Poisson bracket (i.e., is linear in the variables M and p). It is implicit in the 
derivation of this bracket that the integrals are convergent, i.e., the density and momentum 
variables, and the functions of them that appear in the integrals, fall off sufficiently fast 
at large distances. The subscript M indicates that this version employs the momentum 
(as opposed to the velocity) as a dynamical variable. Since several brackets appear below, 
we'll use subscripts to distinguish among them. This bracket, like the versions of the Euler 
equations given above, is unconstrained. It allows for compressiblity, which is, however, 
expressed explicitly only in the Hamiltonian: 



where f = Vx and U represents the internal energy. This bracket and this Hamiltonian 
generate the compressible Euler equations with the pressure given by 

We can apply the bracket ([t]) above to the functionals A4ij,Tjij,i, j = 1,3. Since (see 
equations [s] and [i] above) 

— =x.5,. and — = 0; (10) 
^ = and ^ = x.x„ (11) 
we easily find the following bracket relations: 



{S,„S,,}^^ = (12) 

{Mij,MM}R = SiiMkj - 6jkMii (13) 

{'^ij^Mkijn = + 5jiT.ik. (14) 

These are precisely the relations obtained by Rosensteel by other means ([5], eq. 134); hence 
the index R. Since they are obtained directly from the unconstrained bracket ([T]), they must 
be likewise unconstrained. 
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3. Dirac bracket for the Moment Formulation 

We address here the constraint of incompressibility, which is not incorporated in the 



bracket relations (12)-(14) above. We do this with the aid of the Dirac-bracket formalism. 



We begin this section by defining this and conclude by giving the relations (25), (26) and 



(27) for the Dirac bracket obtained from Rosensteel's bracket. 

3.1. The Dirac Bracket 

Given a bracket {■ , canonical or noncanonical, and an even number 2k of phase- 
space functions {C^}f', one can define a new bracket for which these functions are Casimir 
invariants. This (so-called) Dirac bracket is constructed as follows: 

2k 

{F, G}n = {F, G}-Y, {F, C^} u;-J {C^ G} (15) 

fl,U = l 

where 

u;/^- = {C"^,C'"}; (16) 

it is further assumed that u, an antisymmetric matrix function of the dynamical variables, 
is invertible. The following observations follow directly from this definition: 

1. Each of the functions G^ is a Casimir invariant for the Dirac bracket 

2. Any Casimir invariant of the original bracket {■ , ■} is likewise a Casimir invariant for 
{■,-}d. 



3. {■ , ■}£) is antisymmetric and satisfies the Leibnitz rule (as in equation 20 below). 



It is less obvious but also true that it satisfies the Jacobi identity. This is proved in Appendix 
[Cl below. 

Suppose now that P is a constant of the motion in the dynamics provided by a particular 
Hamiltonian function H under the original bracket, but not a Casimir: {P, H} = but 
{P, G} 7^ for some phase-space function G. Then it is not guaranteed that P will be a 
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constant of the motion in the dynamics provided by H under the modified bracket {■ , -Id! it's 



possible in principle that {P, H^d ^ 0. An example of this is given in Appendix D However, 
some constants of the motion P remain constants of the motion under the modified bracket. 
The following proposition is easily verified: 

Proposition 1. If P = H, the Hamiltonian, or if {P, C^} = for each fi = 1, 2, ... , 2k, 

then P is a constant of the motion in the dynamics provided by H also under the Dirac 
bracket {■ , ■}d- 

In the application of the present paper we find that the constants of the motion are in 



fact unchanged. We have k = 1 and the functions and are given by equation (17) 
below. The only constants of the motion that are not Casimirs of the original bracket are 
the Hamiltionian H and the three components of the angular momentum 

Li = ^ijk-^jk} ^ = 1; 2, 3. 



These commute with and by virtue of the formulas (21) and (22) below. The persis- 
tence of the constants of the motion under the change of bracket follows therefore from the 
Proposition ([T]). 

3.2. A Pair of Constraints 

We choose for the original bracket that of Rosensteel, whose relations are given in equa- 



tions ([12j)-([14j) above. 

As discussed in ^ below, this bracket has a Casimir whose fiuid-dynamical interpretation 
is the magnitude of the circulation vector. In order to construct the incompressible bracket 
we shall augment this algebra by adding two additional Casimirs and expressing the 
constancy of the volume and constancy of the divergence of the velocity field. We may 
express these in the forms 

= In (Det (S)) and = Tr {J^-^M) . (17) 



The explanation for these choices originates in the context of a fluid conflned to an ellipsoidal 
domain and having velocity components that are linear in the cartesian coordinates. Consider 
flrst. The moment tensor S is symmetric and, when transformed to a principal-axis frame 
for the ellipsoid, takes the form 

Q = (m/5)diag (a^, al, Og) (18) 

where Oi, 02, and 03 are the principal- axis lengths and m is the total mas^ Therefore 

det(S) = {m/5f (aiOaaa)^ (19) 

and as deflned above is a constant of the motion as long as the volume (4/8)77010203 is. 
Since for a flgure of uniform density the constancy of the volume implies that of the density, 
the constancy of can be viewed equally as the constancy of the density p. Regarding 
C^, we note that for a fluid having a linear velocity fleld V = L{t)X for some matrix 
L, the divergence of the velocity is the trace of L, which should therefore vanish under the 
assumption of incompressibility. Substituting the expression for V into the moment equation 
(|3|, we flnd that = ^-^M. Therefore 

Tr(L) = Tr(L*) = Tr(S"^A<) 

and the velocity fleld is solenoidal if = 0. 
3.3. Some Useful Formulas 

The calculation of the Dirac-bracket relations and of other related quantities needed below 
requires some preliminary formulas, which we record here. Two useful, general identities for 
matrices A = (Aij) are 

^A^-^ _ .-i.-i , ddetA _ 
— --A,, A, and 



^Thc total mass is conserved; it is a Casimir invariant of the bracket (l7|) (see [T]). 



where Cij is the cofactor of Aij. 

In order to apply the bracket relations to arbitrary functions of functionals, we use the 
derivative propery of brackets: if vi,V2, ■ ■ ■ Vk are functionals and g{v) = g{vi, f 2, . . . , Vk) is 
a real-valued function of them, then for any other functional u 

{u,g{v)} = Y,^{^^^i\- (20) 

1=1 * 

We can now record the following relations for Rosensteel's bracket: 

{S,,,^^}^ = and {S,,,C2}^ = 25,,, (21) 

{M,,C^}^ = -25,, and {M,j,C^] ^ = T.T^Mnj + T.]^Mm (22) 
3.4- Dirac Bracket for Incompressible Ellipsoids 

Since there are only two constraints, the matrix u has only one independent entry, 

= {6,Ai + <5^,S,,) E,/ = 2Tr {^-') (23) 

which implies 



1 
-1 



2Tr(E-i^ 

Thus, the relations for the Rosensteel-Dirac bracket become 



(24) 



{Ei^EkijRD = 0, (25) 
{Mij,Mki}RD = {Mij,Mki}R 

{M.j,^m}b.d = {M.j,^M}R + ,^^0y-_ (27) 

These bracket relations are nonlinear, i.e., the Dirac bracket is no longer of Lie-Poisson 
type. We return later (Q to a verification that they provide a Hamiltonian description of 
the equations governing the motions of the incompressible Riemann ellipsoids. 

9 



4. Dirac Bracket for the Fluid Formulation 



We obtain a different route to tlie Dirac bracket for tlie incompressible, Riemann ellip- 
soids if we first constrain the fluid bracket ([T]) and only subsequently perform the moment 



reduction. For this purpose we carry out the procedure embodied in equation (15) but for 



the original bracket we employ the fluid-dynamical bracket ([T]). We impose the same con- 



straints and as defined in equation (17) and therefore need expressions for the brackets 



{F, C^}m, {F, C'^}m and {C^, C^}a/- For this we need the variational derivatives of and 
with respect to the variables M and p. Straightforward calculations lead to the following: 



5C^ 



and 



and 



-Mki^k^JxiXj. 



(28) 
(29) 



The expressions needed for modifying the bracket are easily obtained with the aid of 



equations (28) and (29): 



JR 



+ 



M 



d 6F 



SF ^3 

R3 OMl 



(30) 



6F 



d X 



(fx, 



(31) 



where A = ^A^S ^; and, from either of the preceding equations. 



{C\C^},, = 2Tr(S-). 



The Dirac-constrained fluid bracket is therefore 



(32) 



1 



2Tr rs 



{{F^ C'}m G}^ - {F, C^}^ {C\ G}J , (33) 



where the index MD denotes the momentum-Dirac bracket, and the index M denotes the 
unconstrained fluid bracket ([7]). We next carry out the moment reduction with the aid of 
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Figure 1: This commuting diagram indicates that the final result is obtained on following either 
path. 



equations (10) and (11). We find 



{S,„C^}^^ = and {S,,,C2},^ = 25,,; 
{M.„ C^},,, = -26., and {M.„ C'}^^ = E-^A^,,. + EJ.'M,.. 



(34) 
(35) 



These are exactly the same expressions as found in the preceding section where the braces 



referred to the finite system of Rosensteel's relations (12), (13) and (14). Since the moment 
reduction of the first term on the right-hand side of equation (33) leads as we have seen to 
Rosensteel's bracket, we arrive at the same constrained, moment bracket via either route, as 
indicated in Figure [ij 

In the next two sections we investigate the structure of Hamilton's equations first using 
Rosensteel's bracket and then using the Dirac bracket based on it. 

5. Dynamical Equations Under Rosensteel's Bracket 

In this section we work out the dynamical equations obtained under Rosensteel's bracket, 
using as the Hamiltonian function the total energy of a Riemann ellipsoid. We shall find 
(see the last sentence of this section) that a key term is missing. 

The symmetric matrix S can be transformed to the diagonal form Q, as in equation 



(18) above. We have implicitly assumed in this description that S is positive-definite: this 
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represents a choice of initial data and, once made, will persist for at least a finite time 
interval. We assign the potential energy appropriate to an ellipsoid with semiaxes ai, 02, 03: 

W = -(1/2) J pV{x)d^x = -(3/10)m^GT, (36) 



where the potential function V is given by equation ([63j) and X by equation (64) of Appendix 
[Aj We have further used equation (22) of chapter 3 of [4j to complete the integration. This 
potential energy is therefore a function only of the squares of the semiaxes, i.e., of the 
eigenvalues of the matrix S. If we write 

E = T*QT, (37) 

we may think of the six independent entries of S as consisting of the three eigenvalues 
together with the three angles needed to specify the rotation matrix T. We may equally 
regard the potential energy as a function of E: 

V(S) = W(Q). (38) 

We use the total energy for the Hamiltonian function: 

H = (l/2)Tr {M'^-^M) + V(S). (39) 

That this function depends only on the moments Ai and S shows that a reduction has been 
achievec|^ . Derivatives of the Hamiltonian are given by the formulas 

,40) 



and 

dH 



dM 



E-^M , (41) 



^This should be compared with Rosensteel's equation (4), where an extra term, proportional to the fluid 
pressure, appears. It is this extra term that leads to the correct dynamical equations under Rosensteel's 
bracket. 
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where indices have been suppressed: to get the ij derivative on the left, one takes the ij 
entry of the matrix on the right. 



We now find, using Rosensteel's bracket relations (12)- (14), the equations of motion 

Mi. = (M'J:-^M).. + ^iM,J:ki}n and (42) 

(43) 



kl 



m + mk 



These represent a dynamical system of dimension eighteen that has a fifteen-dimensional 
invariant manifold expressed by the symmetry of S, and we henceforth restrict considertion 



to this system of dimension fifteen. With the aid of the relation (14) we may rewrite the 
first of these equations as 



dv_ 



(44) 



that 



From the relations = TnQrsTsj (see equation (37) above) and the chain rule, we find 
dW dV dEki dV 



dQij dTjki dQij as 



kl 



-TikTji 



as 



;,j ^^Kl ^^IJ n;6 y / 

where, in the next-to-last term, we have exploited the fact that T and Q may be regarded 



as independent. Equation (44) therefore takes the form 



M = M'J^-'M - 2ST' 



~dQ 



T. 



Next writing M. = T^MT and S = T^QT we obtain the equations in the rotating frame: 



M +[M,Q\=M'Q-^M -2Q^ and Q + [Q, fi] = AT + AT*. 



(45) 



Here f2 = TT* is the antisymmetric angular-velocity matrix and the square bracket is the 



commutator: [A, -B] = AB — BA. At first glance equations (45) do not look like a well- 
determined dynamical system. The variables A/", Q appear on the left-hand side but on the 
right are A/", Q and so this system is well-determined only if f2 is a function of A^ and 
Q. However, the second equation consists of three differential equations for Qn, Q22, Qss 
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and three equations expressing Q = f2(A/', Q). They therefore indeed represent a dynamical 
system of dimension twelve for these variables. The remaining three variables of the original 
fifteen define the rotation matrix T and may be recovered if desired from the equation T = 
QT after the time dependence of Q has been found. Thus we can think of the transformation 
from (I],A4) to {Q,T,Af) as a change of coordinates. 

The equations of motion for the Riemann ellipsoids in their standard form as given in 



equation (62) below likewise represents a twelve-dimensional system. We next bring the 



moment system (45) into this standard form. Recall that the matrix J\f represents the set 



of first moments of the momentum pU where U = Kx is linear in x: 

Mij = j pxiUj £x = J pxiKjiXi (fx = QiiKji. 
Accordingly, we replace the matrix A/" with K through the transformation 

Af = QK\ (46) 



The moment equations (45) take the forms 



K + -nK + Kn = and (47) 

Q = QK' + KQ' + nQ -Qn. (48) 

We focus our attention first on the second of these equations, equation (48). Equation 



(48), which is unchanged under transposition, may be regarded as six equations for the nine 



entries of K. We introduce the matrix 

A = diag(ai, 02, as), (49) 
the matrix of semiaxes, so that Q = {m/5)A'^. The diagonal entries of K are easily found 



to be (for example) Ku = a'l/ai, by virtue of equation (18). Among the off-diagonal entries 



there must be three that are as yet undetermined. If we define a matrix A through the 
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formula 



K = AA-^ + AAA-'^ - fi, 



(50) 



we find tliat equation (48) is satisfied if and only if the matrix A is antisymmetric. This 



prescribes the nine entries of K through the three entries of AA'-^, the three independent 
entries of Q, and the three independent entries of A. This should be compared with 
chapter 4, equation (42), where the same result is arrived at in a different way. 



With the choice (50) for K, we can now express the left-hand side of equation (47) in 



terms of the variables A,A,Q,A. We find 

k + K^ + [K, n] : 



^ + — (AA- riA) + A\-nA + aa^ + n'-A - 2nAA 

dt^ at 



A 



-1 



(51) 



i.e., the left-hand side of equation (47) agrees exactly with that of equation (62) of Appendix 



[a} The right-hand side of equation (47) is diagonal with, for example, the 11 entry 



(52) 



where we have used the definition ( 64 ) of X. 



This gives agreement with equation (62) with the important exception that the pressure 
term is missing. 

6. Dynamical Equations Under the Dirac Bracket 

We address here two aspects of results of the preceding section that are not wholly 
satisfactory. One is the apparent need for a Hamiltonian that is not the total energy as 
usually defined, and the other is that the system obtained is not self-contained but needs 
to be augmented by the further constraints alluded to above regarding the density and the 
divergence. The latter may seem an innocent requirement since such augmentation is needed 
also in the fluid-dynamical derivation as presented in ([3]) or (^); see also the discussion in 
Appendix |A] below. However, the Hamiltonian version, as embodied in the bracket ([T]) above. 
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incorporates not only the law of conservation of momentum but also that of conservation 
of massj^ We should therefore expect the dynamics to be fully described by a Hamiltonian 
description without the need for any augmentation. 



Consider the Dirac bracket {■, - jiji) presented in equations (25), (26), and (27), and again 



employ as the Hamiltonian the total energy (equation 39). The additional terms added by 



the Dirac procedure to the right-hand sides of the bracket relations provide corresponding 
additional terms on the right-hand sides of the dynamical equations. The dynamical equa- 



tions corresponding to equations (42) and (43) therefore become (after a series of tedious 
but straightforward calculations) 

^2 



t = M + M'-2- 



Tr(S-i) 



dQ 



^^^^-^{^-'M+M'J:-') an(B3) 

(54) 



Tr(S-i)' 

where / denotes the unit matrix and is one of the two Casimir invariants of the Dirac 



bracket defined in equation (17) and is therefore a constant of the motion for the preceding 



dynamical system. Since it is proportional to the divergence of the velocity field, it is 
supposed to vanish, and we choose the initial data so that this is so; this simplifies the 
preceding equations. Proceeding as in ^ we obtain from these, with the same definition of 



K as in equation (50) above, the equation 



k + K^ + [K, N] 



+ 



Tr(Q 



{K^ + 2dW/dQ) 



(55) 



This not only has the structure of equation (|62j) but also explicitly provides the expression 
for the pressure that is otherwise obtained by the standard fiuid-dynamical procedure needed 
to maintain the vanishing of the divergence of the velocity field. To see this, observe that 



the term 2pc/p of equation (62) is expressed in terms of the dynamical variables by taking 



^It would also include conservation of energy (or entropy) if we used the full bracket as given in ([!]). 
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the trace of each side of equation (62): 



2pc 
P 



Tr(A-2) = Tr (^if + + Vl] + 2dW/dQ^ = Tr {K^ + 2dW/dQ) . 



Here we have used the identity (51), we have used the formula (52), we have observed that 
Tr{\K^ VL]) =0, and we have set 

T.(/O = |(5:«.M) = |divl7 = 0, 

in accordance with the fluid-dynamical procedure for defining the pressure. This gives for 



the pressure term on the right-hand side of equation (62) the expression 



y = T^^^ Tr {K' + 2dW/dQ) = ^^^^ Tr {K' + 2dW/dQ) Q-\ (56) 
The latter is exactly the extra term provided by the Dirac-bracket formulation and completes 



the verification that the dynamics given by the Hamiltonian (39) under the Dirac bracket is 
exactly that of the Riemann ellipsoids. 

7. Discussion 

Beginning with the Hamiltonian structure of the ideal fluid, we have shown that the in- 
compressible Riemann ellipsoids are governed by Hamiltonian equations in which the Hamil- 
tonian function is the total energy and the constraints of incompressiblity are incorporated 
into a nonlinear bracket via the Dirac formalism. No extraneous constraints are required 
in our formulation. Our results are obtained by introducing a Dirac bracket for the finite- 
dimensional system of moment equations governing the motions of the Riemann ellipsoid, 
and are related in spirit to work of Nguyen and Turski ([?]), who formally introduce a 
Dirac bracket for the purpose of achieving a Hamiltonian formulation of the full, infinite- 
dimensional system of incompressible Euler equations. 

Below we make some additional remarks about constraints. In particular, we show that a 
formulation of Lewis et al. ([8j) for a free boundary liquid, which enforces the incompressibil- 
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ity constraint by requiring divergence free functional derivatives, gives the correct equations 
for a self-gravitating liquid mass. 

7.1. A Bracket for a Free-Boundary Problem 

Lewis et al. ([8J) have proposed the following bracket for a liquid with uniform density 
and a free boundary: 

{FG}-[— (—xu^d^x+f (——-——\d^x (57) 

where uj = curlv and the variations have the following meanings. The functionals F and 
G depend on the velocity field v in the domain D and also on a variable cr determining 
the instantaneous shape of the boundary and defined as follows. The distance Aa is the 
amount that some point x on dD moves normal to itself in the time interval At. Therefore 
(Tt = fi • V. This is the local evolution equation for the motion of the surface normal to 
itself. The variable a is therefore a function of surface coordinates on dD and of time. The 
variational derivative 6F/6\ is clear, but what is less obvious is the requirement that it, like 
the velocity v, be solenoidal: 

div -— = 0. (58) 
ov 

The remaining functional derivative is given by the formula 6F/6(f) = n ■ 6F/6\. It is 
evaluated only on the boundary and is not an independent variation but depends on 6F/6\. 

In their paper, Lewis et al. show how this bracket yields the equations of motion for 
a liquid drop held together by surface tension. We now verify that it does the same if 
surface tension in the Hamiltonian is replaced by self-gravitation. The Hamiltonian is then 
if[v, a] = T[v, a] + W[a], where 

T[v,a]= / {l/2)\Y\^d^x and W[a] = -{1/2) [ Vix)d^x, Vix) = [ 

Jd Jd Jd \^ ~ yl 

The dependence on a arises because the domain D depends on the shape of the boundary. 
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Straightforward calculations show that 



5H 6H ^ 5H , , ,^ , , 



These variations have been made without explicitly imposing the solenoidal constraint (58), 
but note that 6H/ 5v satisfies this constraint anyway by virtue of the solenoidal character of 
V. Therefore 

{F,H}= f ^-{yxlo) d^x+ f h-Y^-^d^x- f n- ^ ((1/2)1 v|2 + \/(x)) rf^x 
Jd civ Jg^ da Jqd c)v 

r SF r 6F 

= / — ■ (-V • Vv + Vr(x)) d^x + / — n-vrf 

Jd civ Jqj;) da 



where we have used the fact that the divergence of 6F/6\ vanishes and a standard vector 
identity. 

On the other hand, 

vj a x + / -^cTj a X. 



Id Jqu 6a 

Hamilton's equations hold if and only if Ft = {F, H} for all functionals F. Comparing 
the expressions for the two quantities we see that we must have = fi • v, expressing the 
free-boundary condition. The equality of the two integrals multiplied by 5F/5\ does not 
guarantee the equality of their coefficients because 6F/6\ is not entirely arbitrary but in 
the Lewis et al. formulation must be constrained by the solenoidal condition: if p{x) is any 
function on D vanishing on dD, jj^ 6F/6\ ■ Vpd^x = 0. Thus the equality of Ft with {F, H} 
implies the correct equation of motion, vj = — v • Vv — Vj5 — W{x), where p is a scalar 
vanishing on dD. 

In principle one should next check whether the moments {I1,A4) effect a reduction with 
the Lewis et al. procedure. Because we know that the Hamiltonian depends only on these 
moments, this amounts to checking that they are closed under the brackets. With the 
definitions of ^ and Q we find for the variational derivatives, on ignoring the solenoidal 
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constraint, 

SMij _ 5Mij _ ^^_n J _ 

ovk OCT ovk oa 

It is seen that SH/Sv satisfies this constraint, but SAi/S\ does not. This can be rectified 
by restricting also the variations 5v to be solenoidal, thereby modifying the expression for 
6Ai/6\ by the addition of a certain gradient. Carrying this out, checking algebraic closure, 
and verifying the equations of motion of the Riemann ellipsoid would require calculations 
of a length and difficulty similar to those already carried out in this paper and we do not 
record these here. 

Fasso and Lewis ([9]) have given an alternative Hamiltonian formulation, not for fiuid 
dynamics, but explicitly for the equations governing the Riemann ellipsoids. 

7.2. The Nature and Number of Incompressibilty Constraints 

The Dirac procedure requires an even number of constraint functions and we have used 
two. It might be surmised that the goal of introducing incompressibility would require only 
one constraint, divf = 0, and that the imposition of a second is an artifice needed in order 
to use the Dirac procedure. This is not so. 



It is easiest to see this in the special context of the Riemann ellipsoids. In equation (62) 
there are two extra parameters, pc and p, that need to be defined in order to make the system 
determinate. One of these is achieved by simply declaring p to be a fixed constant. The 
second is achieved by taking the trace of either side of the equation and setting 

a d 



dt 



diYV = — (^di/a,^ = 0, 



thereby defining function of the velocity field. This definition of pc ensures that 

the preceding equation will hold for all t and therefore that 'Y^di/ai = for all t if this is 
chosen to be true at the initial instant. Our choice of two invariants for the Dirac bracket 
corresponds precisely to these choices. 
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In a more general fluid-dynamical framework in which velocity and density vary with 
position, the imposition of the constraint divf = is not a single constraint, but an infinite 
family of constraints indexed by the position vector x. Once imposed, it implies by virtue 
of mass- conservation equation (|2| that Dp/ Dt = 0, where D/ Dt = d/dt + v ■ W represents 
the convective derivative. This means the initial values of the density are convected by the 
velocity field and necessitates the imposition of a second family of conditions, namely those 
determining the density at the initial instant of time. 

7.3. Invariants 

Notice that the mass m is the zeroth moment of the density distribution and an algebra 
reduction can be constructed for it. It is a Casimir invariant and, as one would expect, so 
is the first moment (the center-of-mass position). By restricting attention to the quadratic 
moments of the density we sit on the symplectic leaf of constant mass and center-of-mass po- 
sition. In the algebra we have constructed, aside from the Casimirs that we have introduced, 
there is one more. 

Rosensteel [3] shows that the magnitude of the Kelvin circulation vector 

= Tr - MM] (60) 

is a Casimir for the algebra {gem (3)) and it remains so for the present algebra]^ That it is 
a Casimir for Rosensteel's unconstrained algebra shows that its validity does not depend on 
incompressiblity. The angular momentum, eij^Mjk is not a Casimir for this algebra, but is 
conserved by the choice of Hamiltonian. 

^This refers to the system in the rotating frame. When the equations of motion are written in the 
inertial frame, it is possible to identify a three-component vector of circulation, each of whose components 
is separately conserved. 
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A. Summary of the Equations Governing Riemann Ellipsoids 

We provide a summary of the basic equation governing the motion of a self-gravitating, 
liquid ellipsoid of spatially uniform density p and semiaxes 01,02,03 with a velocity field 
depending linearly on the cartesian coordinates. A full description is in [1], Chapter 4. 

Relative to a rotating reference frame in which the cartesian coordinates x are aligned 
with the principal axes of the ellipsoid, fiuid motions are allowed that have the form 

u{x) = (^A + AAj A-^x (61) 

where A = diag(ai, 02, 03) and A is an antisymmetric matrix. A and A are in general 
time- dependent, but the full spatial dependence of u is that of linearity in x, as explicitly 
expressed in this equation. The rotation rate of this rotating frame is expressed via a second 
antisymmetric matrix Q, and the dynamical equations governing the time evolution of the 
variables A,Q,A may be written as (cf Chapter 4, equation 57) 

'd'^A d 



+ — {Ak- ^lA) + AX-^IA + Ak' + fiM - 2^1AK 



df^ dt 



A-^ = J-mGA+'^-^A-\ 
2 o 



(62) 

where A = A{A) = diag(^i, ^2j -^s) represents the coefficients in the self-gravitational 
potential 

Vix) = ^mG(l-Y^^Ax^^ , (63) 

which is valid inside the ellipsoid. These coefficients are determined by the semiaxes via the 
formula^ 

f°° du /"°° du I 

1= ——,Ai= where A{u) = J {al + u){al + u){al + u). (64) 

Jo Jo {af + u)A{u) V 

The scalar pc is the pressure at the center x = 0. 



^Thc definitions given here differ by a factor 010203 from those given in (|4]) 
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The system (62 ) consists of twelve first-order equations in the twelve unknowns of A, A, and A 
in which p and Pc appear as parameters. It arises from equation ([T]) only, i.e., from the im- 
position of the law of conservation of momentum only. It must be augmented by further 
information in order to render it determinate. For incompressible flow, two conditions are 
imposed that are consistent with equation ^ of mass conservation: the densitjj^is set equal 
to a constant (which is therefore excluded from the list of variables) and the solenoidal con- 
dition ^dj/aj = is imposed. One can then express pc in terms of the dynamical variables 



A, A, fi, A by taking the trace of each side of equation (62) and putting -^^di/ai = 0; then 
one has twelve equations in twelve unknowns in which the solenoidal condition ^hi/ai = 
is preserved by virtue of the choice of Pc together with the initial datf^ 

B. The Hybrid Coordinate systems 



The transformation of equations (53) and (54) to the equations governing the dynamics 



of ((5,T, A/") was demonstrated in qSl and their equivalence to Riemann's equations of (62) 



with(56) was demonstrated in ||6j Thus Riemann's equations are simply the moment equa- 
tions as we have derived them with velocities and coordinates resolved in a reference frame 
rotating with the body of the ellipsoid. 

In fact, there are four reference frames of interest. The first, with variables 
uses velocities measured in the inertial frame and resolved along axes in the inertial frame 
while the fourth, with variables {Q,N'), supressing the dependence on the rotation matrix, 
measures and resolves velocities in the rotating frame (where Riemann's equations live). 
There is also a second, hybrid frame, with variables (S, A^), where velocities are measured 
in the inertial frame but resolved along axes in the rotating frame and a third, hybrid frame, 
with variables (Q, A/"), where velocities are measured in the rotating frame but resolved along 

^Or alternatively the product aia2a^. 

^Alternatively one can eliminate from the system and achieve a system of ten equations in ten un- 
knowns. 
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axes in the inertial frame. We present the transformation to these frames here. 

For the fourth frame, we showed in qS^that with S = T^QT and J\f = TAiT* the equations 



of motion for (53) and (54), become 



Q=[n,Q]+M + N' and N = [n,N] + N'Q-^N + J' . 



(65) 



where represent pressure and forcing terms which transform in a straightforward fashion. 

In the third frame the velocities are resolved along the inertial frame coordinates but are 
measured along some rotating frame. At the outset there is no need to bias this frame by 
requiring it to be the frame rotating with the body so we can consider an arbitrary angular 
velocity vector uj such that u-^°* = u^"^^* — u; x x. So, defining M.ij = J pXiuf°^ d^x we find 

(66) 



M^j = Mij - J pxiejkiujkxi (fx = Mij - T^u^ij , 
where Qij = eijkUJk- Therefore the dynamics of S and Ai are governed by 



(67) 
(68) 



So far, Cl can be a completely arbitrary, prespecified function of time. The terms on the 



right had side of (68) represent advection, Coriolis, centripetal, Euler and external forces. 



respectively. If we choose a frame to coincide with the body of the ellipsoid, then E must be 
diagonal and, in this manner, Cl is determined. 

The equations for moments completely specified in the rotating reference frame can be 



arrived at by either conjugating (67) and (68) with an orthogonal matrix or by shifting the 



velocity in (65). We shall perform both. Defining A/" = M — QVt and inserting into (65) gives 
easily 



Q = 



(69) 
(70) 
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Alternatively, using Q = TTT^ and Af = TJ\4T^, substituting in (67) and (68), and identi- 
fying Q = TVLT^ gives again (69) and (70). Note, with the above definition of ^2, defining T 



by T = — fiT, results in T* = T. 

C. The Jacobi Identity for General Dirac Brackets 

It is known (cf. [lOj) that a Dirac bracket based on a canonical bracket satisfies the 
Jacobi identity and therefore provides a valid bracket. To our knowldge there is no explicit 
corresponding proof in the literature for the case when the original bracket is more general, 
i.e., not necessarily canonical. We provide that proof here. 

We must show that 

{{F, G}d, H}o + cyclic permutations = (71) 
for all F, G, H and any invertible uj. Therefore 

{{F,G}d.H}d = {{F,G},H}d-{{F,C,}uj-1{C,,G},H}d 
= {{F,G},H}-{{F,G},C'^}uj~J{G^,H} 

C-} {G^ G}, {G\ H] 

where the subscipts on the right hand side have all been dropped in the second line since it 
is unambiuously written in terms of the Lie-Poisson bracket. Upon cyclic permutations, the 
first term will cancel due to the Jacobi identity which holds for the Lie-Poisson bracket, so 
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we can dispose of it immediately. Using the Leibnitz rule, the left hand side of (71 ) becomes 



-{{F, G}, C^} uj-l {C^ H} - {{F, C^}, H} u-l {C^ G} 



-{F, C^} uj-l {{C\ G}, E} - {F, C^}{ uj-l , G} 

+{{F, C"}, CO-', [C^ G} u-l {C^ E} 
G-} u^-l {{G^ G}, < E} 
G^}{ uj-'^ , G^HG^, G} {G^ i/} + c.p.'s. 

The a;~^ term can be pulled out of the bracket in all of the terms by recognizing the relation 

{uj-lF} = -u'lujpliuj^p.F} (72) 



{{G\G^},F}. (73) 



The first three terms and their permutations cancel due to the Jacobi Identity as do the 
second three terms and their permutations. Finally, the last term and its permutations 
cancel amongst themselves due to the Jacobi identity. In this way, it can be shown that the 
Dirac bracket defines a Lie algebra with an even number of Casimirs more than the original 
algebra for any bracket. 

D. Non-Persistence of Invariants 

The Dirac bracket construction ensures that the existence of the Lie-Dirac invariants. 
However, if there exist other dynamical invariants of the unconstrained system, i.e. invariants 
that commute with the Hamiltonian under the unconstrained bracket, canonical or Lie- 
Poisson, then there is no reason that these invariants will remain invariants under the Dirac 
bracket dynamics. Here we give an example where dynamical invariance is lost. 

Consider an A^-body type of system with a Hamiltonian of the form 

N 2 N 2 ^ 

H{P.<D = E| + "^' = E|+E ^(^^ -^.) ' (74) 
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where V{xi — Xj) — V{xj — Xi), and dynamics generated under the canonical Poisson bracket, 

rr X f df dg dg df\ 

y^9^-y[d^^dp, dx^dpj- ^^^^ 



This system conserves the total momentum P — X^^iP?, as is easily shown. 
Now, suppose we constrain away one of the degrees of freedom, by choosing 

C^^xi and C^=pi (76) 

which results in the following Dirac bracket: 

Thus under the constrained dynamics 

dV 

We lose Newton's third law because reaction forces are nulled out by the constraint. 
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